Corner reflectors and Quantum-Non-Demolition Measurements 
in gravitational wave antennae 
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We propose Fabry-Perot cavity with corner reflectors instead of spherical mirrors to reduce the 
contribution of thermoelastic noise in the coating which is relatively large for spherical mirrors and 
which prevents the sensitivity better than Standard Quantum Limit (SQL) from being achieved in 
laser gravitational wave antenna. We demonstrate that thermo-refractive noise in corner reflector 
(CR) is substantially smaller than SQL. We show that the distortion of main mode of cavity with 
CR caused by tilt and displacement of one reflector is smaller than for cavity with spherical mirrors. 
We also consider the distortion caused by small nonperpendicularity of corner facets and by optical 
inhomogeneity of fused silica which is proposed as a material for corner reflectors. 



I. INTRODUCTION 

The existing to-day's multi-layer dielectric coating on 
optical mirrors allows to realize very high resolution exper- 
iments (see e.g. |1J). The reflectivity R in the best optical 
coating has reached the level of (1 - R) ~ 10"^ dH, Q 
(commercially available (1 — R) ~ 10^^), and there are 
many reasons to expect that further improvement of coat- 
ing technologies will permit to obtain the value of (1 — R) ~ 
10"^. With the value of (1 - R) ~ ^0-^ it is possible to 
realize the ring down time Tpp ~ 1 sec in 4 km long Fabry- 
Perot (FP) resonators which are the basic elements in laser 
interferometer gravitational wave antennae (project LIGO 
0,0]). This relatively large value of Tpp permits to have 
relatively small value of the ratio of -^/Xav/Tpp ~ 7 x 1 0^^ 
(if the averaging time Tav — 5 x 10^^ sec). This ratio is 
the limit for the squeezing factor which may be obtained 
if QND procedure of measurement in such FP resonator is 
used 0, Q • Such a procedure will allow to circumvent the 
Standard Quantum Limit (SQL) of sensitivity (see details 
ini). 

Few years ago the role of the thermoelastic noise in 
the bulk of the mirrors was analyzed This analysis 
has shown that if the laser beam spot size on the mir- 
ror surface is sufficiently large then the small value of the 
thermal expansion coefficient asioi — 5 x 10^^ K^^ of 
fused silica will permit to circumvent the SQL sensitivity 
by the factor of ~ 0.1 (if the thermoelastic noise is the only 
source of noises). The consequent analysis of thermoelas- 
tic noise in the coating itself unfortunately predicts that 
the limit of sensitivity will be close to the SQL of sensi- 
tivity lid, lli lllliil- The origin of this obstacle is rela- 
tively big numerical value of thermal expansion coefficient 
ocjaiOs — 5 X 10^^ K^^ of amorphous TaaOs jl3| which is 
used in the best coatings as well as relatively big number 
of layers (usually 20-40) which is necessary to have small 
value of (1 — R). For LIGO project these limitations may 
be illustrated by the following numerical values. The SQL 
sensitivity of detectable amplitude of the perturbation of 
the metric is equal to jl^ 



8h 



2 X 10-24 Hz-^/^ (1.1) 



TUCU^L^ 

where m = 40 kg is mass of test mass, L = 4 km is 
distance between them, cu — 2tc x 100 s^^ is observation 



frequency. Here and below the estimates are calculated 
for numerical parameters listed in Appendix ^ At the 
same time, according to the measurement [T^ , the limit of 
sensitivity of such an antenna only due to the thermoelastic 
noise in the multi-layer TQ2O5 -|- SiOi coating on SiOa 
substrate has to be between [lO| 



;TD coat 



(tu) ~ (0.6^ L4) X 10 



-24 



Hz 



-1/2 



(1.2) 



The goal of this article is to present the analysis of an- 
other version of optical FP cavity where the " contribution" 
of the thermoelastic noise in coating is substantially re- 
duced. The key idea of this version is based on concept 
of corner optical reflector (tri-hedral or two-hedral prism) . 
These types of reflectors were well known among the jewel- 
ers at least from 16-th century (see e.g. autobiography by 
Benvenutto Cellini 0|). In the 70-s of the previous century 
corner reflectors (CR) installed on the Moon allowed to test 
the principle of equivalence for the gravitational defect of 
mass by laser ranging |l(i| |. Here we propose to substitute 
mirrors with finite value of surface curvature (see figH a) 
by 3 facets (figfflb) or 2 facets (see fig^c) corner reflectors 
(CR) manufactured from fused silica. For the same radius 
Rb of laser beam the mass of CR is about the same value 
as cylindrical mirrors (with height about equal to radius 
of cylinder) especially if idle parts of CR is removed. For 
example for 2 facets CR the size of foot surface has to be 
about 45 X 45 cm^ with total test mass about 40 kg (the 
same as planned in advanced LIGO). 

In the proposed scheme the stability of optical mode is 
provided by surfaces lens shaping of each CR foot (as shown 
in figEt^). For reflectors manufactured from fused silica 
the total internal reflection inside the reflectors is possible 
because refraction index nsiOi — ^ -45 is large enough (due 
to Snellius law): nstoi \/2/3 > 1 for 3 facets reflector or 
TLsiOl \/V2 > 1 for 2 facets reflector. 

In sectionn]we consider the modes of ideal cavity (figEJ 
and the distortion of the main mode structure caused by 
small perturbations of different kinds: tilt angle 9 (fig.Ell, 
displacement 5x of one reflector (fig. expose angle e 

(figEJ. 

In section lHll we compare these perturbations for cavities 
with spherical mirrors and with corner reflectors and give 
numerical estimates for the particular case of laser beam 
radius Rb — 6 cm (intensity of beam decreases as e^^ at 
distance Rb from center) which is planned for Advanced 
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FIG. 1: We propose to replace mirrors with finite value of sur- 
face curvature (a) by 3 facets "triprism" type CR (b) or 2 facets 
"roof type CR (c). 



LIGO. Distortions of mode are undesirable in high accu- 
racy spectroscopic measurements because they may pro- 
duce additional noise. For example, in laser gravitational 
wave antenna the light beams from two independently per- 
turbed FP cavities (placed in each arm of Michelson in- 
terferometer) will not produce completely zero field at the 
dark port after the beam splitter. It is equivalent to addi- 
tional noise at the dark port. 

In section Irvl we consider the different sources of optical 
losses of CR and show that they can be at level (1 — R) ~ 
10-5 . 

In cavity with CRs it is necessary to use nevertheless rel- 
atively thin anti-reflective coatings (2-4 layers) on lense 
shape foot. It has to be done to keep the value (1 — R) 
at the level ~ lO^^. Because this coating is substabtialy 
thiner than typical high reflective one (20 - 40 layers) used 
in curved mirrors, thermoelastic noise may be depressed 
by the factor ~ 1 i.e. about one order less than SQL (the 
estimate ()1.2|l is given for 38 layers). The "fee" for use of 
CR is the additional thermo-refractive noise TT] (fluctua- 
tions of temperature produce the fluctuations of refractive 
index) because light beam is traveling inside the corner re- 
flector. However we will show in sectionl^that it is several 



times smaller than SQL for reflectors manufactured from 
fused silica. Important that this thermo-refractive noise 
rapidly decreases with increasing of radius Rb of beam spot 
as - R^;^. 

In cavity with CR the thermoelastic noise in the facets 
remains, but as mentioned above if the reflectors are manu- 
factured from fused silica and the beam spot is large enough 
(see e.g. ^^) then it is possible to circumvent SQL. 



II. THE DISTORTIONS OF MAIN MODE IN FP 
CAVITY WITH CR 

Firstly, we consider FP cavity with two identical per- 
fect corner cube reflectors with three reflecting facets (see 
figCb)- i) the corner angles between the facets are exactly 
equal to 7t/2; ii) the top points of the reflectors are lo- 
cated exactly on common optical axis; iii) the "foots" of 
the reflectors have slight curvatures (shape of a lens sur- 
face as shown in fig. Oi) and they are perpendicular to the 
axis. Then we consider the distortion of mode in this cavity 
caused by different perturbations. 



A. FP resonator with perfect CR 

We can consider that each CR consists of reflector with 
plane foot surface together with spherical lens as shown in 
figl2b- The CR produces mirror transformation, i.e. light 
beam which enters the reflector in point C is transformed 
into the beam leaving the reflector in point C. Using Fres- 
nel integral one can obtain the integral equations for cal- 
culations of eigenmode distribution: 



Go(fl,f2)02(f2)df2 =A(i)i(fi), (2.1) 



Go(fi,?2)Oi(f2)dfi ^X^iiri), (2.2) 
dfi — dxi diji , df2 = dx2 dy2- 

Here functions Oi (fi ) and 02('r2) describe distribution of 
complex field amplitude emitted from imagine flat foot sur- 
faces of reflector 1 (left) and reflector 2 (right) correspond- 
ingly just under lenses in planes AA' and BB'. L is the 
optical path between reflectors (including path inside the 
reflector). Evidently phase fronts coincide with planes AA' 
and BB' so that phases of functions Oi and O2 are con- 
stants. The notation Oi (fi) means the "mirror" transfor- 
mation (produced by 3 facets CR) relative to the optical 
axis^: 

Oilri) = <I)i(-fi), <I)2(fi) = (i)2(-fi) 
The kernel G° is: 



^ ,^ ^ , i i('^^J-T^-Hr(n)-H2(f2)') 

Go(ri ,r2) = - — e V / 

27t 

hi = hi = ^ 



(2.3) 
(2.4) 



^ For 2 facets reflector we have mirror transformation only relatively 
X coordinate: Oi (x,y)) = O] {— x,y) 
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FIG. 3: a). Small tilt of the left CR around its head, b) This 
tilt is equivalent to untilted reflector and displaced lens. 



FIG. 2: a). The perfect alignment of CR assembling cavity, b). 
Each reflector can be regarded as reflector with plane "foot" 
plus a lens. 



Here we use the dimensionless transversal coordinates fi 
and fa (at planes AA' and BB') which can be expressed in 
terms of physical coordinates Ri and R2 as 



b ' 



R2 
b ' 



b = 



k' 



where k is wave vector, hi and ha are additional phase 
shifts produced by spherical lenses at each reflector foot. It 
is easy to see that for spherical lenses the set of eigenmodes 
Oy^"^, OJ^"- and their eigenvalues Amn (below we assume 
Aoo = 1 ) of our FP cavity are described by generalized 
Gauss-Hermite functions: 



1 



A — p 



X exp 

2i(m+n)i4; 



_ =H^ — I X 



(2.5) 
(2.6) 

(2.7) 



-i{m+^/2)^\>--^ 



t]) = arctan 



2r2 



hi(r) =h2(r) = 2r^=2r^ + -3 



Iri = 



(2.8) 
(2.9) 



Here Hralt) is the Hermite polynomial of the order m, Tq 
and Tl are the radii of beam in the waist and at the lens 
correspondingly. It is useful to write down the expressions 
for To, Tl, Tk using g-parameter jl^ (R* — is the radius 
of wave front curvature (in cm) just after the propagation 



of beam through the lens outside of the reflector) : 




sin2\[) = g, cos 2-4) — \/l — g^. 



1-g' 



(2.10) 

(2.11) 
(2.12) 



We are interested in the main mode 0°°(x,'y) of res- 
onator (amplitude distributions of left and right reflectors 
obviously coincide with each other for the main mode: 
(Df(x,y) =(Df (x,ij)). 



B. Distortion due to the Tilt of CR 

Here we consider the main mode 0°°(x,'y) perturbed 
due to tilt misalignment shown in fig. Ot-- We expand 
the perturbed main mode into series over the set of un- 
perturbed modes limiting ourselves to the lowest (dipole) 
approximation: 



«>?°(xi,yi) 



sOOf 



(2.13) 



0^"(x2,y2) + (3f*Ol°(x2,y2). (2.14) 



The tilt of CR around its head through a small angle 
of 6 can be considered as untilted reflector with lens dis- 
placed a small distance 6x ~ 10 perpendicular to optical 
axis (I is the dimensionless distance between the foot and 
head of CR, it is illustrated in fig. ^p) . For this case the 
perturbations of the main mode can be described by dipole 
coefficients defined in (|2. 131 12.1 jjl : 



tilt 



le- 



g(i -g) 



^ V2(1 -g2)3/-4' 

(1-g) [cot (t|;)-i] 
2^/2(^ -g2)i/4 




(2.15) 
(2.16) 

(2.17) 
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FIG. 4: One CR is displaced relative to another one. 

b) 
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FIG. 5: (a). Expose perturbation: the angle between facets of 
corner reflector differs from direct angle by a small value of e. 
It produces the transformation of incident plane wave front (b) 
into a "broken" front of reflected wave (c). 



See details of calculations in Appendix 1X1 

Note that distortion produced by tilt of two facets CR 
around perpendicular axes (angle k on fig. QJ;) can be de- 
scribed by the same formulas as tilt of spherical mirror. 



C. The Distortion due to the Displacement of CR 

One CR can be displaced by a small distance 6x so that 
optical axes of reflectors do not coincide with each other as 
it is shown in fig. 0] For this case the perturbations of the 
main mode can be described by dipole coefficients defined 
by formulae (|2.13l I2.14|) . Denoting the dipole coefficients 

displ J n displ i , ■ 

as a, and (3^ one can obtam: 



(2.18) 



See details of calculations in Appendix IbI 

D. The Distortion of Expose Angle 

Here we consider the case when, for example the left re- 
flector (2-hedral prism shown in fig. ^) has a non-perfect 
perpendicular facets, so that expose angle between them 
differs from 7t/2 by a small angle of e (fig.EK)- Then the 
plane front of incident wave after reflection from the reflec- 
tor is transformed into a broken surface consisting of two 



plane parts declined to the incident wave front by an angle 
of Y = 2e as shown in flg. [S1d,c. 

This statement is also correct for tri-hedral reflectors 
(shown in fig. ^) with the exception of numerical factor: 
if only one facet is declined by angle of e from the normal 
position (and other two facets are non-perturbed) the angle 
Y will be equal to y = 2e-\/2/3. 

Again we can expand the perturbed main mode over the 
set of unperturbed modes of ideal cavity keeping only the 
lowest first-order non- vanishing term of expansion: 



6°° 



of 



, expose -*> 



20 



XI, -yi. 



+ p-P°-(I,f (X2,tj2 



(2.19) 
(2.20) 



(due to the symmetry of this kind perturbation the dipole 
term is null). Calculation gives the following value for aa: 



IyL 



a- 



cxposc I 



4%/2 7Tb^l -g2 

I |^CXpOSC| ^ 



1-g 



Ly 



(2.21) 



(2.22) 



4V27tbg{1 -g2)3/4 
See details in Appendix lO 



III. COMPARISON OF FP CAVITIES WITH CR 
AND WITH SPHERICAL MIRRORS 

Recall that uncontrollable perturbations of the mode 
produce additional noise: in laser interferometer gravita- 
tional antenna the signal at dark port will contain addi- 
tional noise with power proportional to the square of dis- 
tortion coefficients ~ |a|^, |(3|^. In this section we compare 
numerically the distortion of the main mode in traditional 
FP cavity with spherical mirrors (SM cavity) with FP cav- 
ity assembled by CR (CR cavity). 

The distortion of the main mode in SM cavity caused 
by small displacement and tilt of one mirror can be also 
described by coefficients ai of expansion H2.13|l pol . I2H : 



a 



tilt, Sph 
1 



1 



a 



displ, Sph 
1 



V2(1 
(1 



^2)3/4 
,lV4 



V2(1 



g 



1 3/4 



BL 
~b" 

6x 



(3.1) 
(3.2) 



For estimates for both SM and CR cavity we use the pa- 
rameters of cavity proposed for Advanced LIGO j2nl |: 



tl = 6/2.6 ~ 2.3, g = 0.982. 



(3.3) 



These parameters correspond to the radius of laser beam 
Rb at the reflector surface of about Rb — 6 cm. Assuming 
additionally that dimension length I* from foot to top is 
equal to I* = 20 cm, i.e. 

20 cm 



1 = 



7.7, 



and dimension displacement 6x* = b 6x we obtain the fol- 
lowing estimates for SM cavity: 



^tilt, SM^Q_Q13 



10- 



a 



disp, SM 



- 0.0059,^^^ 



6x* 



(3.4) 
(3.5) 
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and for CR cavity: 



IV. THE OPTICAL LOSSES 



^tilt, CR 



1.2 X 10 



-7 



10 



-8 



a 



displ, CR _ 
1 ~ 



0.06 



5x* 



a 



expose 



0.11 



0.1 cm 

(to^) • 



(3.6) 
(3.7) 
(3.8) 



We see that CR cavity is substantially more stable to tilt 
and less stable to displacement than SM cavity. However, 
the total requirements for SM cavity looks more tough if 
one compares the estimates H3.4l and I3.7|l . Indeed to keep 
control of tilt in SM cavity with accuracy 9 ~ 10^^ rad 
(see (jSJ) one has to operate the positioning system with 
accuracy about lb 9 ~ 2 x 1 0^^ cm. But the same level 
of displacement distortion in CR cavity (see (|3.7I) ) one can 
obtain by displacement control with much lower accuracy: 
< 1 mm only (!). 

The requirement for an expose angle in CR cavity looks 
also acceptable: for accuracy of manufacturing e ~ 3x 1 0^^ 
(commercially available prisms have e ~ ±1 x 10^^) and 
hence y — 2ey/2/3 ~ 4.9 x 10^'' we have to take into 
account that three angles between facets (increasing factor 
\/3) in each of two tri-hedral reflector (one more increasing 
factor ^/l) may be independently perturbed (so the total 
increasing factor is equal to V3 x \/2 — \/6) : 



/ exposc\ 



V6afP°^^~0.13 



Optical inhomogeneity is one more source of perturba- 
tion which is specific to CR: the refraction index of fused 
silica (which CR is manufactured from) changes over the 
value of 6n ~ 2 X 10^^ along the length Al ~ 10 cm |F~ 
To estimate negative influence of this effect we can consider 
the model task using the fact that distance scale Al of re- 
fraction index perturbation is about the dimension bl of the 
reflector. Let one half of corner reflector has a perturbed 
refraction index n-|- An(x) which depends on a transversal 
coordinate x only: 



Anfxl = 



An 
An 



— — 6n I 
= 



if 
if 



X > 0, 
X < 



(3.9) 



In this case the beam after reflection from such a reflector 
will have a broken wave front as shown in flg.lsj; with angle 



6n 



(3.10) 



Hence one can estimate the value of perturbation due to 
inhomogeneity using formula (|2.22|) . It is obvious that for 
perturbation with another dependence of An on space co- 
ordinates the formula (|3.1U|) must change this estimate by 
the factor of about unity or slightly larger. So for esti- 
mates we use the 4 times greater value of y than H3.1U|I . 
In that case for two reflector with independently perturbed 
refraction index (factor \/2) we obtain: 



y 



\/2 6n 



2 X 10" 



n 



inhomo 
(X2 



0.011 



(3.11) 



The two last kinds of perturbations (expose angle and 
inhomogeneity) depend only on manufacturing procedure 
and there is a hope they can be decreased due to the im- 
provement of manufacturing culture. 



The loss coefficient for CR cavity must be small — about 
10 ppm. We consider the following sources of losses. 

Fundamental losses on edge are produced by diffraction 
on edge where two facets meet. Qualitatively it can be 
described as two surface waves outside CR (bounded with 
waves inside due to complete internal reflection) meet at 
edge producing diffractional scattering (we acknowledge to 
F.Ya.Fhalili pointed out the existtence of this kind losses). 

To our best knowledge nobody performed rigorous anal- 
ysis of this problem. So we propose the consideration to 
estimate this effect: (a) using the formulas for complex 
coefficient of reflection for plane wave from plane infinite 
boundary between two media for the case of internal re- 
flection (see e.g. '23]) one can construct the solution inside 
CR; (b) using the boundary condition one can obtain the 
fields along outside surface of CR; (c) applying Green's 
formula one can calculate the radiation field in far wave 
zone and total diffractional power. The most vulnerable 
for critics item of this consideration is (a) — applying the 
formulas for infinite boundary to corner configuration. 

We apply this consideration for the case of incident wave 
polarized along edge of CR of "roof type (fig.^^). For this 
particular case (with obvious assumption that magnetic 
permittivity (o. = 1 ) the all field components of constructed 
solution are smooth inside CR and on outside surface CR. 
Our calculation gives the following loss coefflcient: 



(1 -R)d 



0.4A 



0.7 X 10 



-5 



(4.1) 



where A is the optical wavelength, (see details in Appendix 



Note that our consideration of incident wave polarized 
perpendicular to the edge of CR allows to construct smooth 
solution inside CR but this solution on the outside sur- 
face will have break of components of electrical fields at 
the edge. Thus to get a reliable confirmation of the ap- 
proximatetive estimate (14.111 it is necessary either to find a 
rigorous analytical solution of this problem or to perform 
straightforward numerical calculation. 

The losses on non-perfect edge is produced by scatter- 
ing of the plane optical wave on the non-perfect "ridges" 
where two facets meet. The edge of facets intersection with 
uncontrollable width of As < 0.5 |xm will produce optical 
losses which can be roughly estimated as following: 



[1 R] lion— perfect — 



As 



0.8 X lO"'^ 



In other words it satisfies our initial condition to obtain 
T*p ~ Is. 

Optical losses of material. Internal optical losses in pu- 
rified fused silica at the level yioss — 0.5 ppm/cm give 
the loss coefficient about 

(1 - R)opt.ioss ^ yioss2bl ~ 0.8 X 10-5 

Losses in anti reflective coating. As we mentioned in 
Introduction it will be necessary to use anti reflective coat- 
ing on the bottom surfaces of CR. Calculations which we 
omit here shows that to keep the value of (1 — R) at the 
level ~ 10^5 it is sufficient to use 2-4 anti reflective layers 
of coating. 
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V. THERMO-REFRACTIVE NOISE 

The origin of thermo-refractive noise is thermodynamic 
(TD) fluctuations of temperature which produce fluctua- 
tions of phase of Hght travehng inside the CR through de- 
pendence of refractive index n on temperature T: (3 = 
dn/dT 7^ 17]. One can estimate TD temperature fluc- 
tuations using the model of infinite layer with width Ij- 
(0 < z < Ic). We additionally assume that layer is in vac- 
uum and its both surfaces are thermally isolated (thermal 
radiation in accordance with Stefan-Boltzmann law is so 
small that this assumption is quite correct). If light (Gaus- 
sian beam) with radius Rb travels through the layer per- 
pendicular to its surface the fluctuations (p of light phase 
during time t will be defined by TD temperature fluctua- 
tions u averaged over the cylinder ttR^IcI 



cp = klc (3 



(5.1) 



Subscript t means that we are interested in variation of 
temperature during observation time t. 

The total variation of TD temperature fluctuations is 
equal to 



where Tcb is the Boltzmann constant, p is density, and C 
is the specific heat capacity. Then the variation of tem- 
perature over the small time t (adiabatic approximation) 



must be about (u 



(u ) X t/t* where t* 



pCR^ 



/K 



is thermal relaxation time of our cylinder through lateral 
surface (base surfaces of cylinder are thermo isolated), k is 
thermal conductivity. Here we assume that t <C t*. This 
result can be rewritten in form 



pCTtR^lc 




where Vj is thermal diffusive length for the time t. 

Equating (11^)^ — Su(uj) Acu we can obtain the estimate 
for spectral density Su(cu) of averaged temperature putting 
cu ~ Acu ~ 1 /t. 

The accurate expressions for these spectral densities 
Six(cL)) and S(p(cu) of temperature fiuctuations and phase 
fluctuations correspondingly are the following 



Scp(cu) 



4kBT^K 1 



4(32k2lckBT2K 1 



{pC)2 



7TR4tu2 



(5.2) 
(5.3) 



for adiabatic case, i.e. for cu ^ 

One can easy check that our estimate differs from accu- 
rate expression (|5.2|) for Su{cu) only by the factor of about 
unity. In Appendix^ we present derivation of general ex- 
pression for adiabatic and non-adiabatic cases. 

The formulae H5.2I I5.3|l allows to recalculate thermo- 
refractive fluctuations into the fluctuations of dimension- 
less metric h (which usually describes the sensitivity of 
laser gravitational antennae): h, = (p/(kL), where L is cav- 
ity length. It is useful to estimate its spectral density 



Sh(uj) for parameters of laser gravitational antenna (ad- 
vanced LIGO) presented in Appendix IFl and l^. = 10 cm: 



v/ShTcu) ~ 0.5 X 10 



-24 



Hz 



-1/2 



(5.4) 



It is about 4 times smaller than the sensitivity of Stan- 
dard Quantum Limit (II. 1() which is planned to achieve 
in Advanced LIGO 19]. Important that this thermo- 
refractive noise rapidly decreases with increase of beam 
radius: y^StJcuT~ RT ^ . Thus the using so called "mesa- 
shaped" beams |23,|21| (having flat distribution in the cen- 
ter and fall to zero more quickly than Gaussian distribu- 
tion at the edges) with larger radius will allow to decrease 
thermo-refractive noise by several times. For example us- 
ing the 45 X 45 cm^ foot of CR and Rb = 1 cm with mesa 
shape distribution of the intensity of light in the beam one 
may expect the gain of sensitivity for ^ySh(a>) approxi- 
mately one order better than 



VI. CONCLUSION 

The presented analysis of FP cavity with two CR (in- 
stead of mirrors) has to be regarded as an example of cavity 
in which thermoelastic noise in the coating may be substan- 
tially decreased and which permits to circumvent substan- 
tially the SQL of sensitivity and also to have ( 1 — R) ~ 1 0^^ . 

We have shown that CR cavity is considerably more sta- 
ble than cavity with spherical mirrors relative to tilt and 
displacement distortion. The distortion due to expose an- 
gle of CR is not so small but it depends only on manufac- 
turing procedure, and there is a hope to decrease it due to 
improvement of manufacturing culture. 

There does exist another argument in favor of using 
CR in FP resonators. The very recent measurements per- 
formed by LIGO collaborators from University of Glasgow, 
Stanford University, Iowa State University, Syracuse Uni- 
versity and LIGO Lab have shown that multi-layer coating 
also decreases the quality factors of mirrors internal modes 
p^] . This effect may substantially increase the Brown- 
ian component of the noise in the mirror itself and thus 
decrease the sensitivity of LIGO antennae (see details on 
Brownian moise in coating in I25II2II27II2I I2I Isof ) . 

At the same time it is likely that there does exist other 
version of cavity free from thermoelastic noise in coating, 
probably more easy to implement which evidently deserves 
similar in-depth analysis. One of the "candidates" is a 
cavity with unusual reflective coating: each layer of it has 
to have the same small value of thermal expansion as fused 
silica. Unfortunately the technology which may provide it 
is not yet invented. 

In the above analysis we have limited ourselves to the cal- 
culations of the cavity properties itselves and did not dis- 
cuss the coupling of cavity with pumping laser and readout 
system. This analysis has to be done especially because the 
readout system may be an intra-cavity one jq and because 
special attention has to be paid to the possible specific de- 
formation of the main mode distribution in FP cavity with 
CR. One of several possible ways to realize the coupling of 
the mode with pumping source may be based on the exis- 
tence of the evanescent optical field "outside" the surface 
of the facets. In this case it will be evidently necessary to 
use very thin dielectric grating on the surface of the facet. 
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We have not also discuss the polarization characteristics 
of CR. For example it is known that the phase shift of wave 
reflected from plane surface depends on polarization p^ l 
and hence the FP cavity with 2-facets CR will have slightly 
different eigen frequencies for waves with polarization along 
and perpendicular to edge of CR. The additional problem 
to be analyzed is the polarization characteristics of 3-facets 
CR. 

It seems that the CR cavity for Advanced LIGO have to 
be used not for modes with Gaussian distribution of power 
over the cross section but for so called "mesa-shaped" [23, 
I2H mode with flat distribution in the center and fall to zero 
more quickly than Gaussian distribution at the edges. For 
"mesa-shaped" modes the profile of "lenses" (hi and h2 in 
I2.3|l on the foot of the corner reflector must have special 
dependence on radius calculated in [23, 0| . 

It is worth noting that the discussed in this paper fea- 
tures of CR cavity may be useful not only for the grav- 
itational wave antennae but also in other high resolution 
spectroscopic experiments where the low level of optical 
eigenmode fluctuations is important. 



Acknowledgments 

We are very grateful to A. S. Ilinsky and F. Ya. Khalili 
for very fruitful notes and discussions. This work was sup- 
ported by LIGO team from Caltech and in part by NSF 
and Cahech grant PIIY0098715, by Russian Ministry of 
Industry and Science contracts # 40.02.1.1.1.1137 and ^ 
40.700.12.0086, and by Russian Foundation of Basic Re- 
searches grant #03-02-16975-a. 



APPENDIX A: TILT OF CORNER REFLECTION 

The tilt of one reflector through angle (see flgOJi) can 
be considered as displacement of lens on distance 6x ~ 1 9 
(see flgOb) perpendicular to axis. Then the set of integral 
equations (|2.1I |2Jl for eigen mode have the same form with 
replacing kernel Go Gi: 



Go 1 - 



ixi 5x 



(Al) 



The perturbed main modes 0°°(xi,'yi), fl'2°(''2,ll2) we 
flnd as expansion into series over eigen functions of res- 
onator with perfectly positioned reflectors: 

(i)°°(xi,iji) =4)o(yi)X(-iram4>m(xi), (A2) 

m 

<i'2°te,lj2) =4>0(y2)XPm4>m(X2). (A3) 
m 

After substitution these expansions into set (|2.1I O replac- 



ing kernel Go — > Gi wc obtain: 

^A„,(-1)^an,4)^(x2)- 



(A4) 



— e 



ikL 



Go(ri,r2] ^2 I ^ 



X ^(-1)'^am4)m{xi) dxi 

m 

= A^(-1)'^(3^4)m(X2), 



m 

'ixi 5x 



ikL 



Go(fi,r2)x (A5) 

Y |3m4>m{X2) dX2 = 



TR 

Here A is eigen value of perturbed main mode. 
After multiplying equation ljA4|l by cfmo 

(X2) and inte- 
grating we obtain: 



i6xrL 



(A6) 



mo ,0 



V2^ 



After multiplying equation (|A5|) by 4>mo (''i ) integrat- 
ing we obtain: 



AraoCX-Tno — A„^Q_o3mo + Jmo.O i 

i6x Tl 



(A7) 



mo ,0 



V2i 



-^mo-l Pmo-1 a/tOo + 



'>mo + l 



Now we can rewrite equations (|A6IA7|I for different mo 
taking in mind that ao, |3o — 1, 1x1 — Pi = 0(6x), 
^2^(32=0(6x2), ...: 



mo = 0, Ao ao + ^^^^-^"7 Aoai = A Po, 
V2rl 

Aao = AoPo + ^ ^ Al |3i , 
V2ri 



(AS) 
(A9) 



^ao~po^1, A = Ao+0(6x2), (AlO) 



mo = 1 , Al ai — Api 



— iSxTL 

"TIT 



Al , 



A r, — iSxTL ^ 

- Aai + Al |3i ~ ^ T Ao , 



V2i 



>ai 



-iSxtL (a 



f+A^) 



^(3i 



— iSxtL Al 
%/2r2(Ai - Ao)' 



(All) 
(A12) 

(A13) 
(A14) 
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mo =2, A2CX2 - A(32 ^ A2 V2 ai , (A15) 



V2r 



Aa2 + A2(32^ -$^\/2Ai (3i, (A16) 



\/2r 



=^a2 



(A17) 



—^S\r, ^2(^0 c^^ +Ai (3i 

(A18) 



A^-A^ 



Rewriting values o^^ and (3 1 using g-parameter (|2.1UI - I?.12|) 
one can obtain formulas (|2. 151 - 1^. 17|l . 



APPENDIX B: DISPLACEMENT OF CR 

Let the right corner is displaced by value 6x in transver- 
sal direction (see fig. . Then the integral equations for 
perturbed eigen mode is the following 



-"■^^ G] (xi ,yi ,X2,y2) i>^ [^^,V^ ) dx, dyi - 
= Ad2(6x-X2,-'y2), 

Gl (Xi ,111 ,X2,1I2) i>2{~X-2,V2) dX2 d'y2 = 

= Acl)i(-6x-xi,--y,), 

di{-6x-xi,-yi) ~ Oi(-xi,--yi]- 

-6x9x,Oi(-xi -y^], 

02{6x-X2,-'y2) ^ 3'2{-x2,-'y2)+ 

+ 6x9x2 02{-x2 -yi), 
Gi (xi ,yi ,X2,y2) =i Go{xi ,yi ,X2,'y2)x 

X (l -i[5h, +6h2]), 

„ x,6x -X26X 
oh] ~ , 5h2 ' 



(Bl) 
(B2) 
(B3) 
(B4) 

(B5) 

(B6) 
(B7) 

(B8) 



We find perturbed main mode distributions 
^^^(xi ,'yi ), <l'2°(''2,'y2) as expansion into series over 
eigen functions of resonator with perfectly positioned 
reflectors: 

(l)°°(xi,yi) = 4)o{yi)X(-Ta„,4)^{x,) (B9) 

m 

^f[^2,m]^'^o{m]^M^[^2]. (BIO) 



After substitution these expansions into (|B1I - iBTfl we ob- 
tain: 



^ ATT^,o(-l)"^am4>m(X2) - 6 
m 

'i(xi — X2)6x 



ikL 



Go(ri,r2)x 



^^2 ) Y.^-^^"<Xr^'\>r^MAx^ = 

A^{-1)"^|3n,C^n,(x2)+ (Bll) 
m. 

+ A5x^(-1)'-p^(9,,4)^{x2)), 



m 

'i(xi — X2)6x 



kL 



Go(ri,f2)x (B12) 



j ^ |3m4>m(x2) dX2 = 

A Y_ '^m4>m(Xl ) - A5x ^ (9x,4>m(Xl)) , 



After multiplying equation IjBlljl by i^m.oi'^2) and inte- 
grating we obtain: 

Amo.OO^mo.O + Imo.O — A Pmo.O + Jpo.O, (B13) 

2v2rf V 



Iran ,0 



+ ['^mo.O — + 1 ,o] ^/ TUo + 1 amo + 1 j i 

A6x / /mo + 1 _ /mo 

--j:^ I Pmo + l Y 2 Pmo-l Y — 



After multiplying equation (|B12I) by ^rno (''i ) and integrat- 
ing we obtain: 



Amo,o(3mo,0 + Imo.O — Aamo.O + Jpo,0> (B14) 



T' — 
^mn.O — 



iSxTv 



([■^mo.O — Arao-1 ,o] \/molXmo-l + 



j: 



mo ,0 



+ [^mo.O — Am.o + l,o] \/ mo + 1 amo + 1 j . 

-A6x / , ^ \ 

arrio + l V TTlO + 1 - OCmo-1 VTH-O • 



V2 



Here we use the rule: coefficients — and pm = if 

m < 0. 

Substituting Imo,o and Jmo.o into HB13|I and substitut- 
ing I^^ Q and Q into (jB14p we obtain two equations. 
They may be transformed from one to other by replace- 
ment am — > Pm and vice versa pm ^ am- So assuming 
that am = Pm we can solve only one equation: 



'^mo.oamo + _ X 

X (^[Amo.O — Amo-l,o] \ArLoamo-l + 

+ [Amo.O — Amo + 1 ,o] Vmo + 1 
A6x 



(B15) 



=A an 



\/mo + 1 - amo-i \/ttio) , 



We assume that Ao,o = 1 , A = Ao,o + A = 1 -|- A, ao — 1 , 
ai ~ 6x, a2 ~ 6x^, .... Putting mo = in ljB15|) we obtain 
A - 6x^. And putting mo — 1 in IIB15|l we find 



6x / ir 



ai ~ -ao 



+ 



1 



V2rL V2r^ 1 -Ai,o 



(B16) 



Using (|2.1(JI - 12.12|) one can rewrite this formula in form 
(E3HI). 



APPENDIX C: EXPOSE PERTURBATION OF CR 

For this case the equations for eigen modes calculations 
are the following in this section we do not mark by " distri- 
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TABLE L Numerical values of coefficients Fmo.m for low indices 





m = 0Ta = 2m = 4 m = 6 


mo = 


V^Tt 1/\/27t -V(2\/67t) ^/{4^/5n) 



bution function of perturbed mode: 



,ikL 



,ikL 



Go{ri ,r2) f 2(r2) dfa 
AOi (fi ) (j - ibyklxi 

Go{ri,r2)f i(f2)df2 =~ AOilfa). 
We find the solutions as expansion 

m 

^2(^2>'y2) = 4^0(^2) Pm4>m(^2). 

m 

Substituting them into equations (|C1I IC2|) we obtain 



(Cl) 









































^A ^ 



(C2) 
(C3) 

(C4) 



m 

= A ^ am4)m(xi ) (^ - ibyklxi | 

m 

X^n^,0(-1)"''^m*m(X2) =A^(-1)™|3nx4)nx(X2) 



(C5) 



From last equation we obtain 

Pm — ^ 

and substitute it into (jC4p . After multiplying obtained 
equation by (^mo {'X■^ ) and integrating over dxi we get: 

- ibyA^ ktL Y_ '^m T'mo.m, (C6) 
m 

x|4)m(x)4)mo(x) dx (C7) 

We have tabulated coefhcients Fm.mo — result is presented 
in Table n 

Assuming that Ao,o — 1 and ao — 1 , oci , (X2,--- ^ 1 
we see that this system can be divided by two independent 
subsystems: one for odd indices and another one — for 
even indices. Odd indices can be put zero and for even 
indices we have 

iLytL 



FIG. 6: Plane wave traveling and reflecting from dielectric cor- 
ner reflector with angle between facets n/2. The axis z is di- 
rected upward and perpendicular to plane of flgure. Vector E is 
directed along z-axis. 



APPENDIX D: DIFFRACTIONAL LOSSES ON 
EDGE 



Here we write down the calculations to obtain estimate 
(|4.1|) . We consider the monochromatic plane wave polar- 
ized along z-axis traveling and reflecting from dielectric CR 
with angle between facets n/2 (see fig.lSJ. Incident wave: 

Ezinc — to exp ( — icut — iknsin ax — ikncos ay) , 
H=[kE], H = 

Below I drop the multiplier e^^^^. Condition of in- 
ternal reflection is fulfilled on both facets: ncosa > 
1 , n sin a > 1 . We assume that magnetic permittivity 
li. = 1 as in [231 and hence = e (e dielectric permittiv- 
ity). 

Field inside CR. We use Fresnel formulas for light wave 
reflection from plane boundary between dielectric and vac- 
uum using complex reflection coefficient R_l (for the 
case of complete internal reflection): 



Ri(|3) = 



ncos(3 — i-\/TL^sin2(3 — 1 
ncos(3 -I- i\/n^sin^^ — 1 



A^^l 



(X-2 



a4 



bV^' 

iLyktL 



27tb(l 

iLyktL 
2V67tb (1 - ' 



(C8) 
(C9) 
(CIO) 



We see that all coefficients a2, a4, • • • y, i.e. they have 
the same order over y. However the convergence seems 
to take place due to decreasing the coefficients Fmo.o ~ 
1 /ttXq^'^ with rao — > 00. 

These expressions can be rewritten using g-parameter in 

form ^^rniTTT^ . 



where (3 is incident angle. The sum field after two reflec- 
tions from both facets inside dielectric is the following: 

E, = Eo (e-^'^-"-^'^^^' + R^{a)e+'^'^''-'^^y + (Dl) 

+ Ri(7T/2-a)R_L(a)e+^'"''''+^'"«y' 
k=— , kx=krLcosa, ky=knsina (D2) 

The flrst terms in brackets describe the incident wave, the 
second and third terms — the waves reflected from planes 
{y — 0) and (x = 0) correspondingly, the last term de- 
scribes wave double reflected from both planes. 
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FIG. 7: 

We simplify formula I|D1(I and calculate magnetic field: 

= 4Eoe-^'*>'+*«' cos(k,x - 6J cosfk^-y -by], 

(D3) 
V cos^ a — 1 



tan 6y = 



\/ sin^ a — 1 



ncos a 



tan 6y = 



nsm a 



(D4) 



r/ie /le/d on outside surface of CR Now we can write 
down the expressions for fields outside CR in planes x — 
0— e, y = 0— e using boundary conditions — continuity of 
tangent component of E and normal component |i.H. Then 
the expressions for fields are the following: 



Trx=0 
z 

V 



4inEo sin a 6"^'*"+*^ ' cos(6x) sinlk^j-y - 5y ), 
i4nEo cosae^'-'^''+^^' sin(6x) cos(kyi| - 6y ), 
: 4Eo e"^'*" ' cos(k^x - 6^ ] cos(6y ) , 
-4inEo sinae"^'^''+^t'' cos(kxX- 6x) sin(6y), 
-4inEo cosae"^'^''+*«' sin(kxX- 6x) cos(6y) 



We take in mind that plane wave is limited by Gaussian 
multiplier 



b — exp ^ — ([x sin (X — y cos a]^ + 



with small parameter a, or more precisely: a ^ k. 

Radiation field. We can apply diffraction Green for- 
mula to Ez. and calculate it in far wave zone in direction 
characterized by angles cf), and distance |f— f '| (see figtZjl: 



E.ff'l 



G(R1 = 



9nG(R1 



(Ez9nG(R) 

,ikR 



G(R)9nEj df, 



47tR' 



R = If ' - 1^ > - 
k 



-ikG(R) X 



R 



Here f' is radius- vector of observation point, f is radius- 
vector of point on surface of integration, di^ is element of 
integration surface, ft is normal to surface of integration. 



The result of calculations is the following: 
lEo e^^''' e-^^^^+^y'' 



Ez(R) 



ttR 



(D8) 



/7te 



-k^ sin^ 4>/4a^ 



x(l>c + Iy), 



Iy = kb^ 



dx cos(kxX — 6x)e 



ikxcos6 a^x^ sin^ a 



bx (n. COS a sin 6x^1- cos 9 cos 5^ 



(D9) 



ly = kby 



dy cos(ky y - 5y ) e'^^ ^ e-'''^' '=°''' * 
by (n sin a sin 6y — i sin 9 cos 6^ ) 



sin^ a — sin^ I 



(DIG) 



Above we used auxiliary formulas: 



Ic(lc) = 



coskxe dx ^ 7T5(k), if a < k, 



(Dll) 



sinkxe " " dx =^ -, if a < k, 





Power of diffraction losses W can be calculated as fol- 
lowing 



A = 



7t/2 
-7t/2 J 

n 

-n/2 



lEzPc 



-7t/2 



(Ix + Iy)^d0 



"r2 cos4)d4)d9 = |°^ X ^ X A, 
ZTT-^ ka 

(D12) 



2n 

l2 



We have to compare the value W with total power Wq of 
incident wave 

Wo = |^, (1-R), = -^ = ^^xA, (D13) 
2a^ Wo TT^VTtk 

Replacing notations a — > 1 / V2 Rb and calculating numer- 
ically integral A ~ 32.8 for n = 1 .45 (Si O2) and a — n/4 
we finally obtain the estimate H4.1|l . 



APPENDIX E: TD TEMPERATURE 
FLUCTUATIONS IN THERMO-ISOLATED LAYER 

We have thermal conductivity equation for temperature 
u(f, t) in infinite layer with width I (0 < z < I) with fluc- 
tuating force in right part '9] and the following boundary 
conditions: 



= F{f,t), (El) 
= 0, (E2) 

-^^^A6{f-f')6(t-t'), (E3) 



(D5) 


9u 




9t ^ 


(D6) 


9u{f, t) 




9z 


(D7) 




(F(f,t)F{f 



=0, 1 



(pC)^ 

where — K/pC, ke is Boltzmann constant, 6 is Dirac 
delta function. 
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Wc find solution as series: 

dkx dky do) 



u(r,t) 



» (271)3 ^ 
X e^^'t-tKxX-ik^y cosbnZ, 



Ttn 2 1,2 ^2 

bn = — j— , K-j^ — l^-x + l^-y ' 



dxdy dte-^"^*+^'"''''+^'"«^': 



dz ^ cosbnZu(x,'y,z,t), 

L 



We find correlation functions of coefficients Fn(kx, ky , tu) 
Fn.n, = (Fn(l<x,lCy,a))F;;;^ (kx, ,ky, ,a)i)) = 



(pC)2 11^^ 



2-6o, 



X 6(kx-kx,)6(ky -ky,)6(tu-a)i). 

We are interested in temperature u(t,xo,yo), averaged 
over volume V = ttR^I along axis parallel to axis z with 
transversal coordinates Xq and yoi and also its correlation 
function (il(t, 0, 0)il(t + t), xcyo ) with spectral density 
Su(a)): 



u(t,xo,ijo) 



X u(f, t)e 

dz 



dz 



dx di|x 



I 



■°° dxdy 

t) ; 



dkxdkydo) V- „ , , 



X cosbnze^'"'-^'"''''-^'"^^ = 
dkxdkudcu ^i^l 

C 4 X 

. (271)3 

va,t-ik.xo-ik,yo^Q(l^^_^^^)^ (E4) 



X e 



B,n(T) 



;u(t,0,0)u(t + T,xo,ijo) 



IkfiT^K 1 



(pC)2 I 



X e 2- 
kBT2 



dkx dky do) 

k^ e^'"'' 



0)2 + Q4k1 



(pC) 7tR21(1 +2a^r/Rl] 



(E5) 



Sa(uJ1 



dTe^'"''(u(t,0,0)u(t + T,0,0); 



4kBT2K 



k^dki _Ri'<5 
In 



Q4k4 



(E6) 



Making following substitutions 



k^ r2 



2 a-^ pC 



(E7) 



one can express the spectral density using exponential in- 
tegrals: 



S,-,(a)1 = 



ksT^ 
TTpCla^ . 
kRT2 



w2 + £,2 



Ei^(x1 = 



27tpCl 

X (e^^Eii (iw) + e-^'^Ei, (-iw) ) , 
e — 



For particular cases this formula can be simplified: 

4kBT2K 1 



Su(uj)lw<l 



(pC)2l 7tR>2' 



Q i M 4kBT2 r 



PCttR^I r2cu' 



(E8) 
(E9) 



The formulas (|E8I and lE9(l refer to non-adiabatic and adi- 
abatic cases correspondingly. 



APPENDIX F: PARAMETERS 

For our estimates we used the following parameters, ma- 
terial parameters correspond to fused silica. 





= 27tx 100 




A = 1 .064 |xm, 


b 


— 2.3 cm, 


Rb - 


6 cm 


ra 


= 4 X 10"^ 




= 4 X 10^ cm; 


oc 


= 5.5 X 10- 




, K = 1.4xl0 


P 


^22 

cm-^ 


C = 


6.7x10^ i^, 


n 


= 1.45, (3 


dn 
" dT 


= 1.5 • lO"'^ K 



erg 



cm s K ' 
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